Pnicogen-bridged antiferromagnetic superexchange interactions in iron pnictides 
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The first-principles electronic structure calculations made substantial contribution to study of 
high T c iron-pnictide superconductors. By the calculations, LaFeAsO was first predicted to be 
an antiferromagnetic semimetal, and the novel bi-collinear antiferromagnetic order was predicted 
for a-FeTe. Moreover, based on the calculations the underlying mechanism was proposed to be 
Arsenic-bridged antiferromagnetic superexchange interaction between the next-nearest neighbor Fe 
moments. In this article, this physical picture is further presented and discussed in association 
with the elaborate first-principles calculations on LaFePO. The further discussion of origin of the 
magnetism in iron-pnictides and in connection with superconductivity is presented as well. 



I. INTRODUCTION 

Since the discovery of superconductivity in LaFeAsO 
by partial substitution of O with F atoms below 26K[lj, 
intense studies have been devoted to physical properties 
of iron-based pnictides both experimentally and theoret- 
ically. Unlike study of the cuprates, the first-principles 
electronic structure calculations made substantial contri- 
bution to study of the pnictides from the very beginning. 
By the calculations, we first predicted that LaFeAsO is 
an antiferromagnetic semimetal [2j, and then predicted 
that the ground state of a-FeTe is in a novel bi-collinear 
antiferromagnetic order These predictions were con- 
firmed by the later neutron scattering experiments 0, 
Moreover, based on the calculations, we proposed [6[ the 
fluctuating Fe local moments with the As-bridged anti- 
ferromagnetic superexchange interaction to account for 
the magnetism and tetragonal-orthorhombic structural 
distortion in LaFeAsO. 

There are, so far, four types of iron-based compounds 
reported, showing superconductivity after doping or un- 
der high pressures, i.e., 1111-type ifeFeAsO (Re = rare 
earth) [1], 122-type AFe 2 As 2 (A=Ba, Sr, or Ca) [7], Ill- 
type BFeAs (B = alkali metal) [H, and 11-type tetrag- 
onal a-FeSe(Te) [9j. Like cuprates, all these compounds 
have a layered structure, namely they share the same 
structural feature that there exist the robust tetrahe- 
dral layers where Fe atoms are tetragonally coordinated 
by pnicogen or chalcogen atoms and the superconduc- 
tion pairing may happen. A universal finding is that 
all these compounds are in a collinear antiferromagnetic 
order below a tetragonal-orthorhombic structural transi- 
tion temperature 0, OH except for a-FeTe that is in a 
bi-collinear antiferromagnetic order below a tetragonal- 
triclinic structural transition temperature 0, [H, 11 • 

The above finding raises the question on the micro- 
scopic mechanisms underlying the structural transition 
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and antiferromagnetic transition respectively and the re- 
lationship between these two transitions. There are ba- 
sically two contradictive views upon this question. The 
first one [l2| is based on itinerant electron picture, which 
thinks that there are no local moments and the collinear 
antiferromagnetic order is induced by the Fermi surface 
nesting which is also responsible for the structural tran- 
sition due to breaking the four-fold rotational symme- 
try. On the contrary, the second one is based on local 
moment picture. The J1-J2 Heisenberg model was, phe- 
nomenologically [l3[ and from strong electron correlation 
limit [14] respectively, proposed to account for the is- 
sue. Meanwhile, we proposed [6] the fluctuating Fe local 
moments with the As-bridged antiferromagnetic superex- 
change interactions as the driving force upon the two 
transitions, effectively described by the J\- J 2 Heisenberg 
model as well. 

Our proposal embodies the twofold meanings shown by 
the calculations [6[: (1) there are localized magnetic mo- 
ments around Fe ions and embedded in itinerant electrons 
in real space; (2) it is those bands far from the Fermi en- 
ergy rather than the bands nearby the Fermi energy that 
determine the magnetic behavior of pnictides, namely the 
hybridization of Fe with the neighbor As atoms plays a 
substantial role. Here the formation of local moment on 
Fe ion is mainly due to the strong Hund's rule coupling 
on Fe 3<i-orbitals [6j. In this sense, our proposal can be 
considered as the Hund's rule correlation picture. We em- 
phasize again that the Arsenic atoms play a substantial 
role in our physical picture. Subsequently we successfully 
predicted from the calculations [6], based on this Hund's 
rule correlation picture, that the ground state of a-FeTe 
is in a bi-collinear antiferromagnetic order, which was 
later confirmed by the neutron scattering experiments 

Here we would like to address that our physical pic- 
ture works well on all iron-based pnictide superconduc- 
tors, including LaFePO. In this article, we will show how 
our physical picture works on LaFePO through the elabo- 
rate first-principles electronic structure calculations. We 
find that there is P-bridged next-nearest neighbor Fe- 
Fe superexchange antiferromagnetic interaction while the 
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FIG. 1: (Color online) Schematic top view of the FeAs or FeP 
layer in LaFeAsO or LaFePO. The small dashed square is an 
ax a unit cell while the large dashed square is a \[2a x \/2a unit 
cell. The Fe-spins are aligned in a collinear antiferromagnetic 
order. 



nearest neighbor Fe-Fe exchange interaction is very small 
in LaFePO. This results in the collinear antiferromag- 
netic order of Fe moments in the ground state, similar 
to the scenario in LaFeAsO. Moreover, we demonstrate 
that there also exits a small monoclinic lattice distortion 
resulting from this superexchange magnetic interaction, 
again similar to the one in LaFeAsO j6j. 



II. COMPUTATIONAL APPROACH 

In the calculations the plane wave basis method was 
used [15]. We used the generalized gradient approxi- 
mation (GGA) of Perdew-Burke-Ernzerhof [l6| for the 
exchange-correlation potentials. The ultrasoft pseudopo- 
tentials [I?} were used to model the electron-ion interac- 
tions. After the full convergence test, the kinetic energy 
cut-off and the charge density cut-off of the plane wave 
basis were chosen to be 600eV and 4800eV, respectively. 
The Gaussian broadening technique was used and a mesh 
of 16x16x8 k-points were sampled for the Brillouin-zone 
integration. In the calculation, the internal atomic coor- 
dinates within the cell were determined by the energy 
minimization. 

The crystal structure of LaFePO belongs to the tetrag- 
onal structure and the space group is of P4/nmm. In the 
calculations, we adopted two unit cells, namely ax ax c 
crystal unit cell and y/2a x y/2a x c unit cell respectively, 
as shown in Fig. [TJ Here we used the experimental lat- 
tice parameters with a=3.964 A and c=8.512 A. We also 
checked the results with the optimized lattice parameters 
(a=b=3.932 A, c=8.415 A) and there is no significant 
change found. 



III. RESULTS AND ANALYSIS 

In order to derive the exchange interactions between 
the Fe-Fe moments, we designed three different mag- 
netic structures, including the nonmagnetic, the ferro- 



magnetic, and the checkerboard (Neel) antiferromagnetic 
structures calculated in a x a x c crystal unit cell, and 
the collinear antiferromagnetic structure calculated in 
y/2a x y/2a x c unit cell. The calculations show that the 
energy differences among the nonmagnetic, the ferromag- 
netic, and the checkerboard antiferromagnetic states are 
so small that these three states can be considered as de- 
generate, similar to the ones in Ref. [l8[. The calculated 
electronic band structure and the Fermi surface are also 
the same as the ones in the previous calculations [l8|, [l9| . 
From the volumes enclosed by the Fermi surface sheets, 
we derive that both hole and electron carrier density is 
about 2.75 x 10 21 /cra 3 . Thus LaFePO in the nonmag- 
netic phase behaves as a semimetal. However, the calcu- 
lation upon the collinear antiferromagnetic structure in 
LaFePO shows that the ground state of LaFePO is of the 
collinear antiferromagnetic ordering on the Fe moments. 
The energy of the collinear antiferromagnetic ordered 
state is about 33 meV per formula unit cell lower than the 
nonmagnetic one. The magnetic moment around each Fe 
atom is found to be nearly 2.0/i#, smaller than the one 
in the calculation on LaFeAsO [6[. 

These magnetic states can be modeled by the following 
frustrated Heisenberg model with the nearest and next- 
nearest neighbor couplings J\ and J2 



H = J 1 ^2S i 

(ij) 



J 2 X. Si 



Si 



(1) 



whereas (ij) and <C ij ^> denote the summation over the 
nearest and next-nearest neighbors, respectively. From 
the calculated energy data, we find that J\ ~ O.bmeV/S 2 
and J 2 ~ 7.5meV/S 2 . If the spin of each Fe ion 5=1, 
then J\ ^ 0.5meV and J 2 ~ 7.5meV (The detailed calcu- 
lation is referred to Appendix of Ref. [6|) . In comparison 
with LaFeAsO [6], the J\ and J2 are smaller by 7 times 
for LaFePO. 

Fig. [2] (a) and (b) show the calculated charge distri- 
butions on the (001) plane crossing four Fe atoms and 
the plane (110) crossing Fe-P-Fe atoms, respectively. We 
see that there is a strong bonding between Fe and P 
ions. This indicates that the next-nearest neighbor cou- 
pling J2 is induced by the superexchange bridged by P 
3p-orbitals. This superexchange is antiferromagnetic be- 
cause the intermediated state associated with the vir- 
tual hopping bridged by P ions is a spin singlet. The 
charge distribution between the two nearest Fe ions is 
finite but almost none between the next-nearest Fe ions. 
Thus there is a direct exchange interaction between the 
two nearest neighboring Fe spins. Because of the exists of 
the strong Hund's coupling induced direct ferromagnetic 
exchange interaction and the indirect antiferromagnetic 
superexchange interaction bridged by P 3p-orbitals, com- 
bined with the degeneration of different magnetic states 
in a x a x c, the nearest neighbor coupling J\ is nearly 
zero. 

For the collinear antiferromagnetic structure, our cal- 
culations show that there is further small structural dis- 
tortion that shrinks along the spin-parallel alignment (x- 
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FIG. 2: (Color online) Calculated charge density distribution 
of LaFePO in the (001) plane crossing Fe-Fe atoms (a) and in 
the (110) plane crossing Fe-P-Fe atoms (b). 
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FIG. 3: (Color online) Calculated total and orbital-resolved 
partial density of states (spin-up part) of LaFePO in the 
collinear antiferromagnetic state. The Fermi energy is set 
to zero. 



axis in Fig. 1) and expands along the spin-antiparallel 
alignment (y-axis in Fig. 1), causing structural transition 
from tetragonal structure to monoclinic structure. This 
is easily understandable that the direct exchange favors a 
shorter Fe-Fe separation while the superexchange favors a 
larger Fe-P-Fe angle. It turns out that the angle 7 in ab- 
plane (see Fig. [1]) is no longer rectangular and optimized 
to be 90.57° with a small energy gain of about 3meV, 
similar to the finding in LaFeAsO [6j. Such a small dis- 
tortion is found to have a small effect upon the electronic 
structure and the Fe moments. Here we emphasize that 
our calculations show that the driving force upon this 
structural distortion is nothing but the magnetic inter- 
action. More specifically it is the superexchange inter- 
action J2 that breaks the rotational symmetry to induce 
both the structural distortion and the Fe-spin collinear 
antiferromagnetic ordering. 

Fig. [3] shows the total and projected density of states 
of LaFePO in the collinear antiferromagnetic state. The 
most of the states around the Fermi level are of Fe 
3d characters. The corresponding electronic density of 
states is 1.45 states per eV per formula unit cell, which 
is less than the one of nomagnetic state (=5.78 states per 
eV per formula unit cell). The specific heats are evalu- 
ated as 0.84 mJ/(K 2 *mol) and 6.86 mJ/(K 2 *mol) for the 
collinear antiferromagnetic state and nonmagnetic state 
respectively, consistent with the specific heat measure- 
ment [20|. The paramagnetic susceptibility in the non- 
magnetic state is about 1.18 x 10 _9 m 3 jmol. 

Fig. 2] shows the collinear antiferromagnetic electronic 
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FIG. 4: (Color online) (a) The electronic band structure of 
LaFePO in the collinear antiferromagnetic state. The Fermi 
energy is set to zero, (b) The correspondent labels in Bril- 
louin zone, (c) The Fermi surface sheets: a hole-type cylinder 
along rZ and two electron- type pockets between Y and X. 
rX corresponds to the parallel- aligned moment line and YX 
corresponds to the antiparallel-aligned moment line. 
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FIG. 5: (Color online) Calculated total and projected density 
of states at the five Ye-3d orbitals around one of the four Fe 
atoms in the collinear antiferromagnetic state unit cell. The 
Fermi energy is set to zero. 



band structure and the Fermi surface. Unlike the ones in 
the nonmagnetic state [HI , there are now only two bands 
crossing the Fermi level, one small hole cylinder along 
rZ and two small electron cylinders formed between T 
and X. From the volumes enclosed by these Fermi sur- 
face sheets, we find that the hole carrier density is about 
1.34 x 10 20 /cm 3 and the electron carrier density is about 
1.21 x 10 20 /cra 3 . 



IV. ORIGIN OF THE MAGNETISM 

By projecting the density of states onto the five 3d 
orbitals of Fe in the collinear antiferromagnetic state 
of LaFePO (Fig. [5j), we find that the five up-spin or- 
bitals are almost filled and the five down-spin orbitals 
are only partially filled. Moreover, the down-spin elec- 
trons are nearly uniformly distributed in these five 3d 
orbitals. This indicates that the crystal field splitting 
imposed by P atoms is relatively small and the Fe 3d- 
orbitals hybridize strongly with each other. As the Hund 
rule coupling is strong, this would lead to a magnetic 
moment formed around each Fe atom, as found in our 
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calculation. Actually, this is a universal feature for all 
iron pnictides, as we first found in LaFeAsO [6]. The 
polarization of Fe magnetic moments is thus due to the 
Hund's rule coupling. 

In low temperatures, the Fe moments will interact with 
each other to form an antiferromagnetic ordered state. 
These ordered magnetic moments have been observed 
by elastic neutron scattering and other experiments in 
LaFeAsO and other pnictides pi ITlj. However, they are 
not exactly the moments obtained by the density func- 
tional theory (DFT) calculations, as we indicated first 
for LaFeAsO [6|. This is because the DFT calculation is 
done based on a small magnetic unit cell and the low- 
energy quantum spin fluctuations as well as their inter- 
actions with itinerant electrons are frozen by the finite 
excitation gap due to the finite-size effect. Thus the mo- 
ment obtained by the DFT is the bare moment of each Fe 
ion. It should be larger than the ordered moment mea- 
sured by the neutron scattering and other experiments. 
Our calculations show that the bare magnetic moment 
around each Fe atom is about 2/i#. Moreover, the frus- 
tration between the J\ and J2 terms will further suppress 
the antiferromagnetic ordering at the two Fe-sublattices, 
each connected only by the J 2 terms. All of these to- 
gether will reduce strongly the average magnetic ordered 
moment around each Fe measured by experiments, like 
the neutron scattering. 

In high temperatures, there is no net static moment in 
the paramagnetic phase due to the thermal fluctuation, 
but the bare moment of each Fe ion can still be measured 
by a fast local probe like ESR (electron spin resonance). 
Recently the bare moment of Fe has been observed in the 
paramagnetic phase of BaFe2As2 by the ESR measure- 
ment [21]. The value of the moment detected by ESR is 
about 2.2 ~ 2.8 \ib in good agreement with our DFT re- 
sult [g| , which is but significantly larger than the ordered 
moment in the antiferromagnetic phase. 



V. DISCUSSION IN CONNECTION WITH 
SUPERCONDUCTIVITY 

As we see above, the next nearest neighbor antiferro- 
magnetic superexchange interaction affects strongly the 
magnetic structure of the ground state. The magnetic 
fluctuation induced by the antiferromagnetic superex- 
change interactions should be responsible for the super- 
conducting pairing. Upon doping, the long range order- 
ing will be suppressed. However, we believe that the 
remanent antiferromagnetic fluctuation will survive, sim- 
ilar as in cuprate superconductors. We plot the super- 
conducting critical temperatures T c versus J2 for the iron 
pnictide compounds with doping carriers or by apply- 
ing high pressures in Fig. [6j As we see, the maximum 
critical temperatures T c are in proportion to the next- 
nearest neighbor antiferromagnetic superexchange inter- 
action J r 2 . This suggests that there would exist an intrin- 
sic relationship between the superconductivity and the 
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FIG. 6: Various critical superconductivity temperature T c 
versus antiferromagnetic superexchange interaction J2 for 
iron-based pnictides or chalcogenides with do ping or under 
high pressures. Tc are taken from Refs. Q, 043, ^SM- 



Pnicogen-bridged antiferromagnetic superexchange inter- 
actions. 

The above analysis suggests that there are some sim- 
ilarities between pnictide and cuprate superconductors. 
Here we would also like to further address the difference 
between them is that the undoped cuprate is an anti- 
ferromagnetic Mott insulator while the undoped pnictide 
is an antiferromagnetic semimetal. In particular, in un- 
doped pnictide there already coexist both local magnetic 
moments and itinerant electrons. This is similar as in the 
doped cuprate superconductors. In cuprate superconduc- 
tors, the low-energy physics can be effectively described 
by the t — J model, in which the coupling between local 
spins is the antiferromagnetic interaction and the doped 
holes or electrons can hop on the lattice. 

For iron-based pnictides, the low energy spin dynamics 
can be approximately described by an antiferromagnetic 
Heisenberg model (Eq. (pQ) with the nearest and the next- 
nearest neighbor exchange interactions. However, the Fe 
spin (or magnetic moment) is not quantized since the 
electrons constituting the moment can propagate on the 
lattice, like in hole or electron doped high T c supercon- 
ductivity cuprates. Besides these superexchange interac- 
tions, the on-site Hund's rule coupling among different 
Fe 3d orbit als is important. This is because the crystal 
splitting of Fe 3d levels is very small and the spins of 
Fe 3d electrons are polarized mainly by this interaction. 
Thus we believe that the low-energy physical properties 
of these iron-based pnictides can be approximately de- 
scribed by the following effective Hamiltonian 

H = tC if C ia C j(3 ~ J hY1 3i* • Sip 

+ Jl ^ Sia ' Sj/3 + J2 ^ Sia'Sjfi, (2) 

where (ij) and ((ij)) represent the summation over the 
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nearest and the next-nearest neighbors, respectively, a 
and ft are the indices of Fe 3d orbitals. c\ a (c$ a ) is the 
electron creation (annihilation) operator. 

is the spin operator of the a orbital at site i. The to- 
tal spin operator at site i is defined by Si = Si a • In 
Eq. (j2j), Jh is the on-site Hund's coupling among the five 
Fe 3d orbitals. The value of Jh is generally believed to be 
about 1 eV. tf? are the effective hopping integrals that 
can be determined from the electronic band structure in 
the nonmagnetic state [3l|. It has been shown that the 
nearest and the next-nearest neighbor exchange coupling 
constants J\ and J2 in Eq. (j2j) can be calculated from 
the relative energies of the ferromagnetic, square antifer- 
romagnetic, and collinear antiferromagnetic states with 
respect to the non-magnetic state. For the corresponding 
detailed calculations, please refer to the appendix in our 
paper in Ref. 6. Here we have assumed that the contri- 
bution of itinerant electrons to the energy is almost un- 
changed in different magnetically ordered states. Since 
the bare and Jh can be considering independent of 
magnetic structures, the relative energies between differ- 
ent magnetic states are not affected by itinerant elec- 
trons. 



VI. CONCLUSION 

In conclusion, we have presented the first-principles 
calculation of the electronic structure of LaFePO. We 
find that like the situation in LaFeAsO, there are the 
next-nearest neighbor antiferromagnetic superexchange 
interactions bridged by P 3p-orbitals, which are much 
larger than the nearest neighbor ones. This gives rise 
to the collinear antiferromagnetic ordering of Fe spins in 
the ground state as observed in LaFeAsO by the neu- 
tron scattering. However, because of much smaller ex- 
change interactions J2 and J\ in comparison with the 
ones in LaFeAsO, the corresponding transition tempera- 
tures should be very low in LaFePO, very likely less than 
10K scaled according to LaFeAsO. This needs further 
experiment to test and verify. Based on the analysis of 
electronic and magnetic structures, we proposed that the 
low-energy physics of iron pnictides can be effectively de- 
scribed by the t— Jh — J\ — J2 model, defined by Eq. (|2]). 
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